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ABSTRACT 

Massive black hole binaries are naturally predicted in the context of the hierarchical 
model of structure formation. The binaries that manage to lose most of their angular 
momentum can coalesce to form a single remnant. In the last stages of this process, 
the holes undergo an extremely loud phase of gravitational wave emission, possibly 
detectable by current and future probes. The theoretical effort towards obtaining a 
coherent physical picture of the binary path down to coalescence is still underway. 
In this paper, for the first time, we take advantage of observational studies of active 
galactic nuclei evolution to constrain the efficiency of gas-driven binary decay. Under 
conservative assumptions we find that gas accretion toward the nuclear black holes 
can efficiently lead binaries of any mass forming at high redshift (> 2) to coalescence 
within the current time. The observed “downsizing” trend of the accreting black hole 
luminosity function further implies that the gas inflow is sufficient to drive light black 
holes down to coalescence, even if they bind in binaries at lower redshifts, down to 
z ss 0.5 for binaries of ~ 1O 7 M 0 , and z « 0.2 for binaries of ~ 1O 6 M 0 . This has 
strong implications for the detection rates of coalescing black hole binaries of future 
space-based gravitational wave experiments. 

Key words: quasars: supermassive black holes - galaxies: interactions - galaxies: 
nuclei - galaxies: active - black hole physics - gravitational waves 


1 INTRODUCTION 

Massive black hole (M BH) pairs are expecte d to form during 
galaxy mergers (jBegelman. Blandford fe Rees 11 19801 . BBR, 
hereafter). If the nuclei of the two merging galaxies manage 
to su rvive against the tida l forces in play for long enoug h 
(e.g. ICallegari et al. 1 120091 : lYan Wassenhove et al. 1 12014) . 
dynamical friction can efficiently bring the two MBHs to 
the centre of the galactic remnant, forcing them to bind in 
a MBH binary (BHB). 

.(From the binary formation on, dynamical friction 
becomes less and less efficient (BBR), and other dy¬ 
namical processes are needed to further evolve the bi¬ 
nary. In particular, the interaction with single stars and 
wit h nucl ear gas have been thoroughly investigated (see 
iDotti, Sesana fe Decarlill2012l . for a recent review). The sole 
effect of gravitational wave emission forces the two MBHs 
to coalesce within the Hubble time if any physical process 
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manages to shrink the semi-major axes of the BHB down 
to: 


bgw ~ 2 x 10 3 /(e) 


1/4 


A / 4 


(1 + g ) 1 / 2 V 1O 6 M 0 


M 


3/4 


pc, 


( 1 ) 


where M = Mi + M 2 is the total mass of the binary, 
q = M 2 /M 1 is its mass ratio, and /(e) = [1 + (73/24)e 2 + 
(37/9 6)e 4 1 ( 1 — e 2 )~ 7 ^ 2 is a function of the binary eccentricity 
e (IPeters 1 19641 ). The assessment of how effective the various 
processes are to evolve the binary down to bgw is usually 
referred to as the ’last parsec’ problem. 

Attempts to determine the fate of BHBs (whether they 
manage to reach bgw and to coalesce or they remain bound 
in double systems forever) have been made first consid¬ 
ering gas-poor environments, where BHB dynamics is as¬ 
sumed to be driven by three-body interactions with single 
stars. In principle, only stars whose orbits intersect the BHB 
can efficiently interact with it. In an extended stellar sys¬ 
tem, however, only a small fraction of the phase space (the 
so called binary “loss cone”) is populated by such orbits. 
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Stars interacting with the BHB remove energy and angu¬ 
lar momentum from the binary, getting ejected from the 
loss cone. The binary evolution timescale is hence related to 
the rate at which new st ars are fed into the loss cone (e.g. 
iMakino fc Funato 1 1201)4) . Physical mechanisms able to ef¬ 
ficiently refuel the loss cone are required in order for the 
binary to coalesce in less that an Hubble time. Possible 
mechanisms that have been proposed so far are: the pres¬ 
ence o f massive pert urbers (such as giant molecular clouds 
IPerets fc Alexander 1120081 ) . deviations from central symme¬ 
try (e.g. Khan, Just & Merritt 2011; Preto et al. 2011; Gua- 
landris & Merritt 2011, but see also Vasiliev, Antonini & 
Mer ritt 20 1 4) and gravitational potentia l evolving with time 
(e.g. IVasiliev, Antonini fc Merritt Il2014h . 

Similarly, the effect of the interaction between BHBs 
and nuclear gas has been e xplored analytically as well as 
numerically fsee lDotti. Sesana fc Decarli 1120121 . for an up to 
date review). While the full details of the gas/binary inter¬ 
action are still under debate, mainly due to the complexity 
of the system, a clear issue remains to be addressed. Sim¬ 
ilarly to the stellar-driven case, the migration timescale of 
the BHB primarily depends on how much gas is able to spi¬ 
ral toward the MBHs and interact w ith them, instead of, 
e.g. turn into stars le.g. iLodato et al. Il2009tl . This problem 
is remarkably similar to the fueling problem of Active Galac¬ 
tic Nuclei (AGN), i.e. how gas manages to lose most of its 
angular momentum in order to sustain the observed nuclear 
activity. 

Differently from the stellar driven case, however, ob¬ 
servational studies of the AGN population allow to con¬ 
strain the properties of the gas flowing onto a MBH 
(in particular its mass accretion rate). Decades of multi¬ 
wavelength surveys of accreting MBHs have provided a 
relatively robust pictu re of the AG N lu minosi ty func¬ 
tion evolution (see e.g. |H asinger, Miyaii & Schmidt] 2005); 
iHopkins. Richards fc Her cuiist 2007: Buch ner et al. l2015h . 


Coupling such evolution with the observationally deter- 
mined MBH mass function via a continuity equation 
jCavaliere. Morrison fc Wood I Il97ll : ISmall fe Blandford I 
11992h allows to further infer the evolution of the nuclear 
inflow rates as a function of MBH mass and redshift (e.g. 
Merloni & Heinz 2008, Shankar et al. 2013). 

In this work we will assume that the fueling of BHBs 
is consistent with that of single MBHs in the same mass 
and redshift interval. In this way we can estimate the in¬ 
cremental reservoir of angular momentum that a BHB can 
interact with during its cosmic evolution, constraining the 
binary fate, at least in a statistical fashion, directly from 
observations. 


2 MODEL 

2.1 Gas driven BHB dynamics 

To get an order of magnitude estimate of the binary co¬ 
alescence timescale we propose a very simple zeroth-order 
model for the interaction between a BHB and a circum- 
binary accretion disc. We assume that the BHB is sur¬ 
rounded by an axisymmetric, geometrically thin accretion 
disc co-rotating with the BHB. Under this assumption, the 
gas inflow is expected to be halted by the binary, whose grav¬ 
itational torque acts as a dam, at a separation r gap « 2a 


(e.g. lArtvmowicz fc Lubow I [l994h . where a is the binary 
semi-major axis. At this radius the gravitational torque be¬ 
tween the binary and the disc is perfectly balanced by the 
torques that determine the large scale (r a) radial gas in¬ 
flow. Then, we can write the variation of the binary orbital 
angular momentum magnitude (Lbhb) as: 

dLBHB — dLg as — 771 d t \JG A7 T gap (2) 


where M is the total binary mass, and m is the accretion 
rate within the disc. Considering the definition of Lbhb = 
/r%/GMa, where fj, is the binary reduced mass, from eq.[2] 
we derive: 


fi da 
2V2 a 


— 771 dt, 


(3) 


where /r does not evolve in time consistently with the 
assumption that the binary interacts with the disc only 
through the gravitational torques, that stop the gas inflow 
preventing the binary components to accrete. 

As a note of caution, we stress that the circum-binary 
discs could be, in principle , counter-rotatin g with respect to 
the BHBs they orbit dNixon, King fc Pringle |[201l| j. In this 
configuration the gas in teracts with the BHB at ~ a, in¬ 
stead of ~ 2a (|Nixon et al. 1120111 ). and the specific angular 
momentum transfer per unit time is uncertain, depending 
on how strongly the secondary MBH is able to perturb the 
gas inflow. As an example, if all the gas passing through 
a bounds to the secondary MBH, the binary angular mo¬ 
mentum diminishes by two times the angular momentum 
carried by the gas. In this scenario, the BHB evolves on the 
same timescale regardless of the BHB-d isc relative orienta¬ 
tion. Moreover, iR.odig fc Sesanal (j2014h demonstrated that 
BHBs embedded in self-gravitating retrograde discs may sec¬ 
ularly tilt their orbital plane toward a coplanar prograde 
equilibrium configuration. For these reasons we will focus 
our investigation on the prograde case in the following. 

As a first order of magnitude estimate of the coalescence 
timescale we can make the simplifying (although quite com¬ 
mon) assumption of an Eddington limited accretion event 
and integrate eq.QT| Further assuming a fixed radiative effi¬ 
ciency (e = 0.075, see section 2.2) we obtain 


AfBHB 


/Oi\ fiec 2 

Va c / 2V2L Edd 


10 ' 


(i + g) : 


■In 


(I) 


yr (4) 


where a; and a c define the binary separation range where 
the MBH-gas interaction drives the binary orbital decay. 

Eq.lH shows that, in order to coalesce, the binary has 
to interact with an amount of matter of the order of its 
reduced mass, with only a weak dependence on the exact 
ratio between the initial and final separation. A conservative 
estimate for this ratio can be obtained setting m equal to 
the radius at which the two MBHs bind in a binary 


di ~ GM/2ol ~ 0.5 ' PG (5) 

where we have assumed the M — relation (iGultekin et al. I 
2009). The final separation can be conservatively estimated 
as a c ~ 6 x 10 _5 (M/10 6 Mq) 3 ^ 4 pc, in order for the BHB 
to coalesce due to the emission of gravitational waves in 
~ 10 4 yr. Under these assumptions a,i/a c oc M and 
ln(ai/a c ) < 9 for a binary mass > 10 6 Mg. 
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2.2 Gas inflow onto BHBs: observational 
constraints 


We can now, for the first time, try to relax any a priori 
assumption on the accretion rate in the circum-binary disc 
(such as the Eddington limit used in eq. [4j), assuming an 
observational driven prescription for the evolution of m as 
a function of MBH and redshift. In particular, we adopt the 
average accretion rates obtained assuming that the MBH 
evolution is governed by a continuity equation, where the 
MBH mass function at any given time can be used to pre¬ 
dict that at any other time, provided the distribution of 
accretion rates as a function of black hole mass is known. 
The continuity equation can be written as: 

^ (' ip(m,t)(M(m.,t ))) = 0 (6) 

where m = Log M ( M is the black hole mass in solar 
units), ip(m,t) is the MBH mass function at time t, and 
( M(m,t )) is the average accretion rate of a MBH of mass 
M at time t. The average accretion rate can be defined 
through a “fueling” function, F(rh,m,t), describing the dis¬ 
tribution of accretion rates for objects of mass M at time t: 
( M(M,z )) = f MF(m,m, z) dm. Such a fueling function is 
not known a priori, and observational determinations thereof 
have been able so far to probe robustly only the extremes 
of the overall population. However, the fueling function can 
be derived by inverting the integral equation that relates 
the luminosity function (rf>) of the AGN population with its 
mass function. Indeed we can write: 




F(£ — m, t) dm 


(7) 


where we have called t = Log Lbol and £ = Log (ec 2 ), with t 
the radiative efficiency. This is assumed to be constant and 
its average v alue can b e estimated by means of the Soltan 
argument (ISoltan Ill982l ~). which relates the mass density of 
remnants MBH in the local Universe with the integrated 
amount of accreted gas during the AGN phases, as identified 
by the luminosity function. 

Gilfanov & Merloni (2014) reviewed the most recent 
assessments of the Soltan argument. Adopting as a starting 
point the bolometric AGN luminosity function of Hopkins 
et al. (2007), the estimate of the (mass-weighted) average 
radiative efficiency, (e), can be expressed as: 

« 0.075 Ko(l -ii- Cct + Ciost)]" 1 (8) 

where G) = Pbh,z=o/4.2 x lO 5 M0Mpc~ 3 is the local 
z = 0) MBH m ass density in units of 4.2x 10 5 A/gMpc -3 
Marconi et ah 2004); G is the mass density of black holes 
at the highest redshift probed by the bolometric luminosity 
function, z k, 6, in units of the local one, and encapsulates 
our uncertainty on the process of BH formation and seed¬ 
ing in proto-galactic nuclei; £ct is the fraction of the MBH 
mass density (relative to the local one) grown in heavily ob¬ 
scured, Compton Thick AGN; finally, Goat is the fraction of 
BH mass contained in “wandering” objects, that have been 
ejected from a galaxy nucleus, for example, in the after- 
math of a merging event because of the anisotropy in the 
emiss ion of gravitational waves (e.g. iLousto fe Zlochowerl 
120131. and references therein). More recent estimates of the 


fraction of MBH mass densit y accumulated i n heavily ob¬ 
scured, Compton-Thick AGN jBuchner et al. 1(201 oi l suggest 
that ^ct ~ 0.35. Neglecting G i n ec f HI (he- assuming a 
negligibly small seed BH mass density), the average radia¬ 
tive efficiency will vary approximately between 0.075 and 
0.1 for 0 < Gost < 0.3. Therefore in the following we will 
use the results obtained by performing a numerical inver¬ 
sion of eq. [3 based on a minimization scheme that used the 
Hopkins et al. (2007) AGN bolometric luminosity function 
as a constraint, and assuming a fixed radiative efficiency in 
the range 0.075 < e < 0.1. The average Eddington ratios 
/Edd (bolometric luminosity normalized to the Eddington 
limit) and accretion rates obtained in this way are shown as 
a function of redshift in the left and right panels of figure [T] 
respectively. We note that increasing the radiative efficiency 
value implies a decrease in the average accretion rates, espe¬ 
cially for higher MBH masses at higher redshifts where the 
MBH evolution is relatively more important. This is consis¬ 
tent with the adopted calculation scheme where the AGN 
luminosity function is assumed as a constraint to derive the 
accretion rates estimates. 


3 RESULTS 

The observational constraints on the average value of rh 
(function of M and z) discussed in section 2.2 allow us 
to numerically integrate eq. [3] to determine the migration 
timescale AfBHB for any BHB. We can further translate 
AfBHB into an estimate of the minimum redshift zbhb at 
which a BHB of mass M must form in order to coalesce 
within a given redshift z co al- The results of the numerical 
integration are shown in figure [3 for a binary mass ratio 
q = 0.1 (lower panel) and q = 0.3 (upper panel). No ’last 
parsec’ problem seems to exist for binaries of total mass 
M < 10' Mgand q < 0.3 formed at zbhb > 0.5. More 
massive binaries (M > 10 s ) M@ do not coalesce within the 
present time if formed at 2bhb < 1-2, and the extreme cases 
of M > 10 9 Mq coalesce within z = 0 only if formed at 
zbhb iG 2. Binaries forming at higher redshift coalesce in 
shorter times, since the average accretion rates increase with 
2: within the redshift interval considered in this analysis. For 
example, all the q < 0.3 binaries forming at 2bhb ~ 2.5 co¬ 
alesce within 2 coa i = 2. 

As a note of caution, we stress that the assumption 
that all the gas inflow is stopped by the binary is over¬ 
simplifying. As a matter of fact, numerical 2-D and 3-D 
simulations (independently of the exact treatment of grav¬ 
ity or hydrodynamics) demonstrated that the deviations 
from axisymmetry close to the binary, driven by its grav- 


r, a „ (lHavasaki. Mineshiee & Ho 2008; Cuadra et al. 

2009; 

Rbdig et al. 

201ll: Sesana et al. 20121; Shi et al. 

2012; 

Noble et al. 

2012 

; D’Orazio, Haiman & MacFadven 

2013; 

Farris et al. 

2014 

). To put firm upper limits to the MBH 


migration timescale, we assume that only a fraction / = 0.4 
of the gas inflow interacts dynamically with the binary, 
while the remaining 60% of the gas fails to strongly inter¬ 
act gravitationally with the binary, and falls onto one of the 
MBHs unimpeded. A simple timescale estimate for the bi¬ 
nary shrinkage, under the assumption of MBHs accreting at 
a fixed fraction of the Eddington limit, can be obtained re- 
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Figure 1. Average Eddington ratios (left panel) and mass accretion rates (right panel) of MBHs as function of z. Black, red, green 
and blue colors refer to MBH masses of 10 6 , 10 7 , 10 8 , and 10 9 Mq, respectively. The shaded areas show the range of values comprised 
between the two limiting cases considered for the radiative efficiency (see discussion in the text) corresponding to e = 0.075 and e = 0.1. 



Figure 2. Minimum redshift zbhb a t which two MBHs must bind 
in a binary of total mass Mbhb in order to coalesce within three 
different redshifts z coa \ = 0 (black lines), 1 (blue lines) and 2 (red 
lines). The upper and lower panels correspond to a mass ratio 
q — 0.3 and q = 0.1, respectively. As in figure [I] the shaded areas 
correspond to the range of values obtained considering a fixed 
radiative efficiency comprised between e = 0.075 (left curve) and 
e = 0.1 (right curve). Thin dotted lines mark the three values 
assumed for z coa \, and are shown to facilitate the reader. The 
case represented in the figure is when all the gas inflowing toward 
the galaxy nucleus interacts with the secondary MBH. 


Figure 3. Same as figure [2] either assuming that all the gas in¬ 
teracts with the secondary MBH (shading with inclined lines) or 
that only a fraction of 0.4 of the inflowing gas mass interacts with 
M 2 (shading with horizontal lines). 


placing Z/Edd with / Ledd in eq. [4] Reducing the fraction of 
interacting matter increases the migration timescale, hence 
increasing the minimum zbhb required for the BHB to coa¬ 
lesce within a given z coa 1 as shown in figure 0 

As expected, decreasing / the evolution of every binary 
slows down, but the general trends discussed while com¬ 
menting the f = 1 cases remain valid. BHBs with total mass 
M ~ 10 7 M©at z > 2 are particularly affected, because of 


© 0000 RAS, MNRAS 000,[I][6] 






































Observational constraints on MBH binary pairing 5 


the redshift at which their typical Eddington ratio peaks 
(see figure [I]). Still, these BHBs manage to coalesce between 
Zcoal ~ 1 — 2, as well as their more massive counterparts. 


cone of the binaries at rates significantly higher than those 
expected in spherica l system s j, possibly leading BHBs to a 
fast coalescence fsee lVasiliev 12014. for a recent discussion). 


4 CONCLUSIONS 


ACKNOWLEDGMENTS 


We estimated the gas driven orbital decay of BHBs from the 
instant at which they bind in a binary down to their final co¬ 
alescence. For the first time we propose an observationally 
driven approach, that has the advantage of not being af¬ 
fected neither by any assumption on the (largely unknown) 
feeding process driving the accretion, nor by the fraction of 
the gas inflow that turns into stars at large scales before 
interacting with the BHB. 

Our investigation proves that 1) high redshift BHBs of 
any mass coalesce on a very short time; 2) Low mass BHBs 
(M < 1O 7 M 0 ) formed at low redshift manage to merge any¬ 
way within z = 0, since their accretion history peaks at 
lower redshifts. These findings are particularly relevant since 
the coalescence of low mass BHBs is one of the sources of 
gravitational waves detectable by future space based grav- 
itation al wave interferometer s, such as the mission concept 
eLISA llAmaro Seoane Ifeoi.lh . 

We have worked under very conservative assumptions: 

- We assumed that binaries in the late stages of galaxy 
mergers are fueled as much as MBHs in comparable isolated 
galaxies, without assuming any merger driven boost in 
the accretion. The merger process itself is considered, 
however, an efficient reshuffler of the gas angular momen¬ 
tum at galactic scales, driving efficiently gas inflows all 


vations (e.g. 

Kennicutt & Keel 

1984 

Keel et al. 

1985; 

Alonso et al. 

20071: Koss et al. 

20 111: 

Ellison et al. 

2011; 

Silverman et al. 201ll: SatvaDal et al. 

2014|l as well 

as by 


a wealth of nume rical work s performed on differen t kind 
of merger s fe.g. iDi Matte o, Springel & Hernquist I 20051 


Johanssonl Burkert fe Naab 20 091: Hopkins fe Quataert 


2010 ; ICallegarietaL 201inVan Wassenhove et al. 1 120121 ; 

Capelo et al. I l2014h . 


- We have assumed that all the gas accretion onto MBHs is 
radiatively efficient, and that the mass accretion rate at few 
gravitational radii (R g , where basically all the luminosity 
is emitted) is equal to the that at thousands of R g , where 
the gas interacts with the secondary MBH. We stress that 
a significant fraction of the accretion flow, however, could 
be ejected in the form of fa st outflows, a s often found in 
numerical simulations (e.g. iProga I 120031 : iNaravan et al. I 

l2012ft . 


Under our conservative assumptions, gas driven migra¬ 
tion of high mass (M > 10 8 Mq) BHBs formed at low red¬ 
shift could be inefficient. Such binaries are of particular in- 
terest, b eing the onl y ones observable through pulsar timing 
(IHobbs et al. 1 I2010I . and references therein). The morpho¬ 
logical and dynamical characteristic of their hosts suggest, 
however, that interactions with stars could play a significant 
role in the binaries shrinking. The hos ts of very massive 
MBHs often show tria xial profiles (e.g. iFaber et al. IH997I : 
iKormendv et airil2009l . and references therein). The lack of 
spherical and axial symmetry in the potential of the hosts 
allows the single stars to modify substantially their angu¬ 
lar momentum components. Stars can hence re-fill the loss 


We acknowledge the anonimous Referee, Alberto Sesana and 
Eugene Vasiliev for useful comments and fruitful discussions. 
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